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In order to fully characterize the noise associated with electron transport, with its severe consequences 
for solid-state quantum information systems, the theory of full counting statistics has been developed. It 
accounts for correlation effects associated with the statistics and effects of entanglement, but it remains 
a non-trivial task to account for interaction effects. In this article we present two examples: we describe 
electron transport through quantum dots with strong charging effects beyond perturbation theory in the 
tunneling, and we analyze current fluctuations in a diffusive interacting conductor. 
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1 Introduction 

Solid-state quantum information systems based on electronic spin or charge degrees of freedom offer a 
number of intrinsic advantages and drawbacks. Among the former are the fast operation times, the possi- 
bility to scale the systems to large size, and the relative ease to integrate them into electronic control and 
read-out circuits. Probably the most serious disadvantage is the fact that, in general, solid state devices 
suffer strongly from noise due to internal and external degrees of freedom as well as material-specific fluc- 
tuations. They lead to relaxation and decoherence processes. Hence, one of the major tasks in the field is 
the understanding and control of noise and decoherence. In this article we will concentrate on the analysis 
of fluctuations which arise due to the discrete nature of the electron charge. They lead to what is denoted 
as shot noise. Its analysis reveals information about electron correlations and entanglement (TJ. 

In many circumstances the fluctuations are Gaussian distributed and fully characterized by their power 
spectrum. In order to get further information in the general case, e.g. about correlations and entanglement, 
the theory of full counting statistics (FCS) 1 2 1 of electrons has been developed. It concentrates on the 
probability distribution for the number of electrons transferred through the conductor during a given period 
of time. It yields not only the variance, but all higher moments of the charge transfer as well, and thus also 
contains information about rare large fluctuations. 

The FCS has its historical roots in quantum optics, where the counting statistics of photons has been used 
to characterize the coherence of photon sources |3 1. Photons detected by the photo-counter are correlated 
in time, reflecting the Bose statistics of the particles involved. For electronic currents the Fermi statistics 
is the relevant one, but the first attempts to derive the FCS of electrons |4 1 revealed some fundamental 
interpretation problems, related to subtleties of the quantum measurement process. We are interested in 
the probability that the outcome of a measurement of the charge is q. According to text-book definitions 
of projective measurement this quantity can be expressed by (n\S(q — Q)\n), where Q is the operator for 
the transmitted charge, and \n) denotes the quantum state. It is tempting to relate the transferred charge to 



Copyright line will be provided by the publisher 



4 



D. A. Bagrets, Y. Utsumi, D. S. Golubev, and G. Schon: Full Counting Statistics of Interacting Electrons 



the current operator, Q — f*° j(r)dr. However, in quantum transport problems one has to pay attention 
to the fact that the electric current j(r) is an operator, which, in general, does not commute at different 
times. This property led to severe interpretation problems within the original work. They were resolved 
in the later work of Levitov and Lesovik [2 1, which invokes explicitly an extra degree of freedom, namely 
the detector degree of freedom. The paradigm of projective measurement is then applied to this detector 
degree of freedom. 

In the meantime the theory of FCS in mesoscopic transport has developed into a mature field; some 
achievement are summarized in Refs. (0151. However, the experimental analysis of the FCS remains a 
challenge. First measurements of the third cumulant of charge transfer through a tunnel junction have been 
reported in Refs. |6|[7] an d, very recently, the FCS of a semiconductor quantum dot (QD) has been inves- 
tigated by a real-time detection of single-electron tunneling via a quantum point contact 1 8 1. Furthermore, 
threshold type of measurements of the FCS using an array of over-damped Josephson junctions has been 
elaborated theoretically in Ref. [9|. 

In this paper we will review our recent results on the FCS of interacting electrons in a QD and low- 
dimensional diffusive conductors. The QDs are basic constituents of most solid-state quantum information 
systems. For example, superconducting single Cooper-pair boxes 1 10 1 and, similarly, a double-dot system 
formed in a semiconductor 2DEG [ 1 1 1 have been shown to operate as charge qubits. A metallic QD or 
single-electron transistor (SET) can serve as an electro-meter to measure the quantum states of a charge 
qubit 11 21 . Since all these devices are based on the charge measurements, a thorough understanding of the 
fluctuating properties of charge becomes crucial for progress in this field. 

Recently further links became apparent between the FCS of electron transport and the field of solid-state 
quantum information processing. One of these is related to the use of electron entangled states for these 
purposes. Most of the work on entanglement has been performed in optical systems with photons 1131 . 
cavity QED systems 1 14 1 and ion traps 1 15 1. By now several ideas have been put forward how to generate, 
manipulate and detect electronic entangled states 1 16 1. It turns out that in solid state systems entanglement 
is rather common, the nontrivial task remaining its control and detection. For mesoscopic conductors, the 
prototype scheme of such detection was discussed in Ref. 1171 . It has been shown that the presence of 
spatially separated pairs of entangled electrons, created by some entangler, can be revealed by using a 
beam splitter and by measuring the correlations of the current fluctuations in the leads. If the electrons are 
injected in an entangled state, bunching and anti-bunching of the cross-correlations of current fluctuations 
should be found, depending on whether the state is a spin singlet or triplet. In Ref. |18| the FCS of 
entangled electrons has been analyzed in detail. The FCS depends not only on the scattering properties of 
the conductor but also on the correlations among the electrons that compose the incident beam. In Ref. 1191 
the Clauser-Horne inequality test for the FCS in the multi-terminal structures has been proposed in order 
to detect the entanglement in the source flux of electrons. 

A second link is the intrinsic relation between FCS and detector properties of a quantum point contact. 
QPCs were suggested as charge detectors in Ref. |20| and have been studied experimentally in Ref. |21 1. 
Recently they have been used as detectors for the state of quantum-dot qubits I22l l23ll24l . The operating 
principle of the QPC detector relies on the dependence of the electron current / through the QPC on the 
state of the two-level system. In Ref. 1 25 1 the detector properties of the QPC have been calculated beyond 
linear-response for arbitrary energy-dependent transparency and coupling. This is the case of interest since 
for maximum detector sensitivity typical measurements are done in the regime of high QPC transparency, 
D ~ 1/2, and for coupling that is not weak 1221 . It was found that both the back-action dephasing rate T 
and the measurement rate W are determined by the electron FCS. 

A further motivation to study FCS arises from the need to understand the effect of interaction on elec- 
tron transport in disordered low-dimensional conductors. Disorder enhancement of Coulomb interaction, 
together with quantum coherence effects strongly influence the transport properties of these systems. The 
FCS analysis of this long-standing problem provides a deeper insight into the question. Typically the 
Coulomb interaction leads to a suppression of the conductance of mesoscopic samples at low temperatures 
and bias voltages. It has been demonstrated 1 26l 1271 l28l that the strength of this suppression in various 
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types of mesoscopic conductors is related to their noise properties. Frequently we find the simple rule: 
the higher the shot noise the stronger is the Coulomb suppression of the conductance. The reason is that 
both shot noise and Coulomb corrections to the transport current are manifestations of the discreet nature 
of the electric charge. Beyond that, it was shown that the Coulomb correction to the shot noise scales with 
the third moment of the current in the absence of interactions |29|. Furthermore, it was demonstrated by 
renormalization group studies of the FCS of short coherent conductors [ 30 1 and of quantum dots 1311 in 
the presence of Coulomb interaction that the interaction correction to the n~ th moment of the current is 
determined by n + 1-th moment evaluated in the absence of interaction. Further developing these ideas 
we show in the present paper that Coulomb interaction may substantially enhance the probability of large 
current fluctuations in low dimensions, leading to the appearance of long correlated 'trains' in the trans- 
ferred charge. This effect is most pronounced when the system size matches the dephasing length due to 
Coulomb interaction. Such coincidence is not accidental and comes from the presence of the soft diffusive 
modes in the system which strongly renormalize the bare interaction. 

The structure of this article is as follows: In the next section we introduce some basic definitions and 
concepts of the FCS in mesoscopic transport. We discuss the paradigm of quantum measurement by using 
a spin 1/2 as galvanometer and consider some simple illustrative examples. In the main part we concen- 
trate on our own contributions to the field, discussing the effects of Coulomb interaction onto the shot 
noise and FCS in interacting quantum dot systems (Section 3) and in low-dimensional diffusive interacting 
conductors (Section 4). 



2 Concepts of FCS 

We start this section by introducing some definitions and general formulae of the FCS approach to meso- 
scopic transport. The central quantity is the probability distribution, P to (N), for N electrons to be trans- 
ferred through the conductor during a time interval to- The detection time to is assumed to be much larger 
than the inverse current frequency e //, which ensures that on average N ^> 1. This probability distribution 
P to (N) is related to the cumulant generating function (CGF), T(x), via a discrete Fourier transform 



e 

N=— 



Pt (N)e m \ (1) 



The auxiliary variable \ is usually called "counting field". From the CGF one readily obtains the "cumu- 
lants" (irreducible moments) 



Ck = {(N k ))=-(-i) k ^(x)\ x=0 . (2) 
The first four of the irreducible moments, defined by 



Ci = N = J2NP t0 (N), C 2 = (N~N)i, C a = (N-N)Z, 



N 



C A = (N — N) 4 — 3C%, (3) 

denote the mean, variance, asymmetry ("skewness") and kurtosis ("sharpness"), respectively. They charac- 
terize the peak position, width of the distribution and further details of the shape of the distribution Pt (N), 
as illustrated in Fig.Q 

In order to provide a quantum mechanical definition of the CGF of electrons we will follow the approach 
proposed by Levitov and Lesovik 1 2 1 . The key step is to include the measurement device in the description. 
As a gedanken scheme a spin- 1/2 system is used as a galvanometer for the charge detection. This spin 
is placed near the conductor and coupled magnetically to the electric current. Let the electron system be 
described by the Hamiltonian 7i(q, p). We further assume that the spin-1/2 generates a vector potential 
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Fig. 1 The distribution of the number of transmitted electrons N. The mean G\, the variance C2, the skewness C3 
and the kurtosis C4 characterize the peak position, the width, the asymmetry and the sharpness of the distribution, 
respectively. 



a(r) of the form a(r) = ^X^ f( v )- Here the function /(r) smoothly interpolates between and 1 in the 
vicinity of the cross-section at which the current is measured, and \ is an arbitrary coupling constant so 
far. It will be shown below that it plays a role of "counting field". If one further restricts the coupling 
of the current to the ^-component of the spin then the total Hamiltonian of the system takes the form 
Ha = W(q, p - &a z ). 

In the semiclassical approximation, when the variation of V/(r) on the scale of the Fermi wave length 
\p is weak, it is possible to linearize the electron spectrum at energies near the Fermi surface. Thus one 
arrives at the Hamiltonian H a = H(q, p) + Hint* where 

1 r+°° v 
Hint = —a z d 3 r a(r)j(r) =-±a z I S - (4) 

e J-00 2 e 

Here j(r) is the current density and Is — J d 3 r j(r) V/(r) the total current across a surface S. On the 
quasi-classical level Eq. (|4} shows that a spin linearly coupled to the measured current Is (t) will precess 
with a rate proportional to the current. If the coupling is turned on at time t = and switched off at t 
the precession angle 9 = X Jo" Is(t)dt/e of the spin around the z-axis is proportional to the transferred 
charge through the conductor. In this way the spin- 1/2 turns into analog galvanometer. 

To proceed with the fully quantum mechanical description let us consider the evolution of the spin 
density matrix ps(t). We assume that initially the density matrix of the whole system factorizes, p = 
Pe <8 />s(0), with p e being the initial density matrix of electrons. Then the time evolution of p(t) is given 
by 

p(t) = Tr e (e-^ t pe^ t ), (5) 

where Tr e denotes the trace over electron states. Since, by construction, the evolution operator e~ %u,,t is 
diagonal in the basis of a z , the spin density matrix takes the form 



PS (t ) 



Pn (0) Z(x)p n (0) 
Z(-x)p^(0) Pu (0) 



Z(x) = Tre(e- in * t Pee* n -x t ), (6) 



where the Hamiltonian TL X = W(q, p) Is acts on the electron degrees of freedom only. It becomes 

— ( ' 

clear now that the non-diagonal elements of the density matrix (|6) contain the information about the distri- 
bution of precession angles of the spin during time to. To make it explicit we use the transformation rule 
of the spin- 1/2 density matrix corresponding to a rotation around the z-axis by the angle 9, 



Kg(p) 



Ptt « ' PU 



e PIT PU 



(7) 



Copyright line will be provided by the publisher 



pop header will be provided by the publisher 



Fig. 2 Keldysh contour C, used for evaluation of the cumulant generating function (CGF) . 

One now can identify Z(x) with the CGF introduced in Eq.Q, i.e one sets Z(x) = e^^^, and the spin 
density matrix ps(to) can be represented as a superposition of the form 

+00 

ps(t )= Pt n (N)n e = Nx (p), (8) 

JV=— 00 

where P to (N) has a meaning of the probability to observe the precession at angle 9 = N\- For a classical 
spin a precession angle 9 = \ corresponds to a current pulse carrying an elementary electron charge, e = 

to 

Is(t)dt. Using the correspondence principle we conclude that the quantity P to (N) can be interpreted 



as the probability of transfer the multiple charge Ne. This consideration suggest the use of Z{x) as the 
microscopical quantum mechanical definition for the CGF. 

Using a cyclic permutation under the trace, Tr e , in Eq. one can represent Z(x) in the form of the 
Keldysh partition function 

Z{ X ) = e~ nx) = (T K ™p\-i I dtU x {t)\\. (9) 



Here the time integration is performed along the Keldysh contour C, as shown in Fig.|2] and Tk denotes the 
time ordering operator along the path C. The average (...) is performed with the non-equilibrium electron 
density matrix p e . The interaction part of the Hamiltonian H x (t) reads Wi n t(t) = ^jxW-Tss where the 
"counting field" x(^ ± ) = iX i s asymmetric on the upper and lower branches of the Keldysh contour. 

The definition l|9} for the CGF can be generalized to obtain the full frequency dependence of the current 
correlators of arbitrary order |32|. Consider a mesoscopic conductor as shown in Fig. [5] coupled to two 
leads such that lead 1 is grounded while lead 2 is biased with voltage V(t), We assume that the current I(t) 
is measured in the lead 2. This set-up is described by the interaction Hamiltonian Tt- m t(t) = ip(t)I(t)/e 
with the phases y ± (t) = /_ eU(r)dr ± defined on the lower/upper branches of the contour C, 

and Eq. (|9j yields the generating functional Z[(p + (t), <p~{t)\ of the current fluctuations in the lead 2. In 
analogy with the definition the higher-order derivatives of the functional yield the n-point irreducible 
correlation function of currents 

e^Cnih, ...,t n ) = He)"^ . . . \nZ[ V +{t)^-{t)]\ x=0 . (10) 

For illustration we consider an Ohmic resistor with resistance R at temperature T. Its CGF is quadratic 
in (p^ 1 (t) and reads 

r R [<P + (t),ip-(t)} = ^-?§- [ dh [ dh a{h - faM*iM*2), (ID 



2tt R 



c Jc 



where Rq = 2-Kh/e is the quantum resistance and 

nT 2 

a(t) = ^ . (12) 

V ; 2sinh 2 (7rfT) 

If the times ii,2 lie on different branches of the contour C the kernel a(t) is regularized by the shift into 
the complex plane, — > t ± iO; otherwise it is understood as principal value. The quadratic form of 
reflects the Gaussian nature of current fluctuations in an Ohmic resistor with its well-known properties. 
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R 




Fig. 3 Voltage biased mesoscopic conductor of resistance R. 



The simplest example of a CGF for non-Gaussian processes is provided by a tunnel junction. In this 
case one has 1331 



-4?r^ / dti [ dt 2 a(h 
Rt Jc Jc 



ti) sin 2 



(13) 



where Rt is a tunnel resistance. For a constant applied voltage eV and stationary "counting field" x the 
corresponding CGF reduces to 



F T (eV, x) = -to [r+ (e 4X - l) + T_ (e~* x - l)] , r± = ±- 



1 



eV 



2 R T 1 - e ±eV 7 T ' 



(14) 



This result represents the CGF of a bidirectional Poissonian process with rates T± corresponding to un- 
corrected tunneling processes of the charge through the junction. At zero temperature and positive bias 
voltage the second term of Eq. J 1 4b disappears. Then, performing the inverse Fourier transformation we 
obtain a simple Poisson distribution, corresponding to uncorrected charge transfer 



P(N) = 



« dx 



exp(J r T ( e y,x) - iNx) 



N\ 



N = t r 



(15) 



As a further important example we consider a quantum point contact (QPC), shown in Fig.|4] It can 
be used as a quantum detector for the state of a quantum dot charge qubit. Its operating principle is based 
on the property that - due to the electrostatic coupling between the dot and the QPC - the scattering matrix 



Si 



t* 



and thereby the current / of the QPC depend on the state \j) of the qubit. For a given 



realization S of the scattering matrix the FCS has been calculated by Levitov et al. | 2!|, with the result 

de 



-to 



2tt 



In det 



1 + / s+(- x )S(x)-i 



(16) 



where S(x) — ( , i X /2 5 * J , and / = diag(/L, is the diagonal density matrix of the leads. 
The explicit evaluation of Eq. (II 61 yields 



HX) = -to J ^ In [1 + f L (l - f R )D(e*x - 1) + f R (l - f L )D{e~^ - 1)] 



(17) 



where D = \t\ 2 is a transmission coefficient. The physical interpretation of this result is that electrons 
can be transmitted either forward or backward with probabilities pr^l = — fn)D and pn-R — 
/fl(l — Jl)D, respectively, with the occupation factors accounting for the Pauli principle. 

Recently, it has been realized that the quantum detector properties of the QPC are intrinsically related 
to its FCS |25|. The two basic quantities to be considered are the measurement-induced dephasing time 
of the qubit, and the time needed for the acquisition of information about the qubit's state. Keeping in 
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Fig. 4 The principal scheme of the qubits readout using the quantum point contact (QPC). Due to electrostatic 
coupling the current / driven by the external voltage V is controlled by the state \ j) of the qubit. 



mind the gedanken spin- 1/2 measurement scheme for charge detection, described in the beginning of this 
section, we realize that in the quantum measurement process by a QPC the role of gedanken galvanometer 
is played by the qubit (or more generally by a many-level system). Then the analog of Eq. © describes 
the dephasing rate T of the qubit's density matrix p(t) due to its interaction with the electron current in the 



QT>C, p jk (t) = p jk (0) 



giHi-tg—iHjt 



Pjk(0)e 



. Here the average is taken over a stationary state of 



the QPC and H, and H k are electron Hamiltonians describing the propagation of electrons with scattering 
matrices Sj and S k . In the long-time limit, t > %/eV, the decay is exponential with rate 



- jk 



de , , 
— In act 

2tt 



i + fists, 



1) 



(18) 



The analogy with the expression (1161 is striking. For a qubit at low temperature, T <C eV, one obtains 1 25 1 



ri2 = -^ln|<i^+r 1 r*| 



(19) 



The second aspect of the measurement by the QPC is the rate of information acquisition. The infor- 
mation about the state \j) of the qubit is encoded in the probability distribution (N) for N electrons 
to be transferred via the QPC, given that the qubit is in the state The mean value of this distribution, 



Nj, and its width y (AN 2 ) grow like t and t 1 ! 2 . This time dependence implies that only after a certain 

time, which we denote as the measurement time r mcas , the two peaks, corresponding to different states | j) 
and | k) emerge from a broadened distribution. Quantitatively this time can be defined by considering the 

statistical overlap of two distributions, M jk (t) = J2n {N)Pjf^ ( N ) ■ For lon S times, t > h/eV, 
the decay should be exponential, Mjk(t) cx exp{— Wjkt}, with Wjk — l/^mcas being the measurement 
rate. As it was shown in Ref. 1 25 1 it can be expressed in terms of the CGF as 



w ik = ttt min x [Fjix) + Fk{-X)\ 



(20) 



In the case of quantum-limited detection the rates Wjk and Tjk coincide, while generally Wjk < T^, 
meaning that dephasing occurs faster than the information gain. One can show that a QPC can be operated 
as quantum limited detector with the rates 



eV 



where R j(k) = 1 - D j(k ) 



(DjD k )V 2 + {R.Ru) 1 ' 2 
are reflection coefficients. 



(21) 
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(b) 



-1/2 _ , +1/2 
Qo/e 



Fig. 5 (a) The equivalent circuit of a metallic quantum dot. (b) The Coulomb oscillations. 



The FCS of tunnel junction and quantum point contact represent the two simplest generic examples with 
non-Gaussian current fluctuations. In both cases the effects of electron-electron interaction were neglected. 
This is justified provided the conductance of the system is not too small, G > e 2 / 27r?i, and the length (size) 
of the system does not exceed the inelastic mean free path. Accounting for interaction effects on the FCS 
is no trivial task. This is the case, in particular, when electrons have an internal dynamics within the 
conductor, so that in the derivation of F(x) me internal degrees of freedom have to be integrated out. Two 
examples of these calculations are considered in the following sections. In Section 3 we discuss the effects 
of Coulomb interaction onto the shot noise and FCS in interacting quantum dot systems, and in Section 4 
we address the statistic of current fluctuations in the low-dimensional diffusive interacting conductors. 



3 Full Counting Statistics in interacting Quantum Dots 

In this section we review the FCS of a single electron transistor (SET), which is shown schematically in 
Fig. |5](a). It consists of a QD with strong local Coulomb interaction, which is coupled via two tunnel 
junctions (left and right) with tunneling resistances and _Rr and low capacitances Cl and Cr to two 
electrodes (source and drain). It is further coupled capacitively via Cq to a gate electrode which allows - by 
an applied gate voltage which determines the 'gate charge' Qq — CqVq - to control the number of elec- 
trons inside the QD. In the systems of interest the total capacitance Cs = Cl+Cr+Cg is low and, hence, 
the single-electron charging energy Ec — e 2 /2Cs typically in the range of IK or above. If Ec exceeds the 
temperature and applied source-drain bias voltage, | eV\, the electron transfer through the QD is suppressed 
("Coulomb blockade"). However, the Coulomb barrier can be tuned by the gate charge. For example, the 
energy difference between the charge zero state and the one with one excess electron in the QD depends 
on Qq as Ao = Ec(l — 2Qq/b). The voltage drops across the tunnel junctions are /i r = n r eV, where 
k l/r = = ' = Cr/l(Cl + Cr )~ 1 and the subscript r stands for either L or R junction. If these voltages satisfy 
the condition /iR < Ao < /iL electrons can tunnel sequentially through the island. This mechanism leads 
to the typical oscillating behavior of the conductance as a function of the gate charge illustrated in Fig.|5](b). 

The first steps toward a theory of FCS in quantum dots with account of Coulomb interaction has been 
performed in Ref. 1 34 1 where the FCS of charge pumping in the limit of high transmission of the contacts 
was considered. Further progress was made by one of the authors and Nazarov |35|, who derived the 
FCS in the frame of a Master equation. This approach is valid in the weak tunneling regime, where the 
parameter g = Rq/Rt, i.e., the ratio between the quantum resistance Rq and the effective (parallel) 
resistance, i?^ 1 = R^ 1 + , is much smaller than unity, g -C 1. Applied to a quantum dot in the 
vicinity of the first conductance peak, the CGF is found to differ from a simple Poissonian distribution. 
Rather it reads 



^Xx)=t,T ^-\ D( X ) = l + l^(e--l) + ^i(e---l). (22) 
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Here T = Til + Tir + IYi + Tri, and the rates of electron tunneling into/out of the island through the 
junction r are given by Fermi's golden rule, 

For special cases the CGF can be simplified: (i) Close to the threshold of the Coulomb blockade regime, 
e.g., for A « /il, and low temperatures T <C \eV\ the tunneling process through junction L becomes the 
bottleneck since Tli oc — Ao) is much smaller than Tir. In this case the CGF reduces to a Poissonian 
form, 

~i r L i(e lx -l) . (24) 

(ii) For a symmetric SET, i?L = Rr and Cl = Cr , at the conductance peak, Ao — 0, one finds for T = 
and eV > 

JT« „ 2 N (e^ 2 - 1), e 7V/i = V/2(R L + R R ) . (25) 

The extra factor 1/2 in the exponent leads to a sub-Poissonian value of the Fano factor, i.e., ratio between 
dc power spectrum and average current Su/2e(I) « 1/2, indicating that tunneling processes through the 
two junctions are correlated. The distribution function in this case becomes P(N) =X}jv n =o ^ p (-^l) ■ 
Pp(N-r) <5Ar,(Ar L +Ar R )/2' where the distributions of N-^ and TVr transmitted electrons through the L and 
R junctions have a Poissonian form Pp(N) = N N e~ N /Nl . Both are constrained as indicated by the 
Kronecker 8. 

The Master equation approach captures the basic physics of the strong Coulomb correlations inside the 
QD, but it neglects non-Markovian effects, which become important for strongly conducting QDs, i.e., 
if the dimensionless conductance g is no longer small. This includes quantum fluctuations of the charge 
due to co-tunneling, i.e., simultaneous tunneling of two electrons through two junctions. This process 
dominates in the Coulomb blockade regime, i.e. far away from the conductance peaks. Recently Braggio 
et al. 1 36 1 considered these effects in second order perturbation theory in g in extension of the theory 1 35 1 
using well established real-time diagrammatic techniques 8371 1381 . 

The CGF of a quantum dot in the limit of very strong tunneling, g 1 , has also been considered 1 31 • In 
this limit the Coulomb blockade almost disappears. Its weak precursor is caused by quantum fluctuations of 
the phase, which is the variable canonically conjugated to the island charge. The small negative correction 
to the conductance is logarithmic: g w g — 2 In(f2/T), where Sl = 1/{RtCy) is the inverse RC time |39|. 

Several further articles dealt with different setups. In Ref. |40| bosonization techniques were used to 
find the FCS of an open quantum dot coupled to reservoirs by single-channel point contacts in the presence 
of a strong in-plane magnetic field. Similarly, the CGF for the generalized two-channel Kondo model, 
which models a QD in the Kondo regime, has been derived |41 1. In both cases the authors succeeded to 
fully account for Coulomb correlations, but the results are limited to a very special, exactly solvable case. 
Despite this work the understanding of the effects of quantum fluctuations on the FCS of interacting QDs is 
far from complete. In what follows we evaluate CGF for the regime of intermediate strength conductance. 



3.1 FCS OF A SET FOR INTERMEDIATE STRENGTH CONDUCTANCE 

Here we consider a quantum dot single-electron transistor in the intermediate strength tunneling regime, 
where (introducing for convenience a new dimensionless conductance parameter) ao = g/(2ir) 2 < 1, We 
assume that the inverse RC time is still smaller than the characteristic charging energy, f2 <C Ec, which 
ensures that the charge-state levels are well resolved. In the vicinity of the conductance peak, precisely for 
\Ao\/Ec 1, it is sufficient to restrict the attention to only two charge states of the quantum dot with 
charges differing by e. The Hamiltonian can then be mapped onto the 'multi-channel anisotropic Kondo 
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model' 1421 . Introducing a spin- 1/2 operator a, which acts on the charge states, we have 

H = ^ '^£rka] kn a rkn + ^-a z + ^ X! { T ^\ kn a Yk > n a + + H.c). (26) 

r— L,R,I kn r— L,Rfcfc'n 

Here a\ kn creates an electron with wave vector k and channel index n (including spin) in the left or right 
electrode or island (r=L,R,I). The tunneling matrix elements T r are assumed to be independent of k and n. 
The junction conductances are 1/R T = 2ire 2 N cri \T r \ 2 pip r , with iV c h being the number of channels and p T 
the electron DOS. We assume that energy and spin relaxation times are fast so that electrons are distributed 
according to a Fermi distribution function f(u>) = l/[exp(ui/T) + 1] both in the island and in the leads. 

In the intermediate conductance regime the main consequence of the quantum fluctuations of the charge 
is the renormalization of system parameters, specifically of the charging energy and the conductance. 
A perturbative two-loop renormalization group analysis for ^> 1 predicts a renormalization of the 
conductance, oq — » zo^o, ar, d of the charging energy, Ao — > zoAq, to depend logarithmically on the low 
energy cut-off A = max{T, A } ll42l 

Zo= l + 2a ln(£c/A)' (2?) 



Such a conductance renormalization has been confirmed by experiments 1 4-3 1, where the observed height 
of the conductance peaks has been suppressed as 1/ In T. 

The logarithmic renormalization is typical for Kondo problem, in which one encounters logarithmic 
divergences in perturbation theory. Likewise, a perturbative treatment of quantum fluctuations in a quantum 
dot leads to logarithmic divergences. Handling these divergences remains a nontrivial task, especially in 
non-equilibrium transport problem. Schoeller and Schon 1 37 1 have formulated a real-time diagrammatic 
approach to this problem. Summing up a certain class of infinite order diagrams, they managed to remove 
the divergences and recover the renormalization factor J27> . Furthermore, they derived non-linear current- 
voltage characteristics including low-bias Kondo anomalies. Recently the second cumulant of the current, 
i.e. the noise, has been evaluated in lowest |44| and second-order perturbation theory |45 1. However, apart 
from the second-order analysis by Braggio et al. 1 36 1, the FCS of a quantum dot in the moderate tunneling 
regime has not been yet analyzed, in particular in situations where the finite-order perturbation theory fails 
and infinite order diagrams need to be included. Motivated by that, we addressed this problem in Ref. 1 46 1, 
where the FCS of a SET has been evaluated with the use of Majorana fermion representation 1 47 48 49 1 . 
This formulation enabled us to apply Wick's theorem and consequently the standard Schwinger-Keldysh 
approach 1501 I51ll52l . Since practical calculations are rather technical, we will first summarize our main 
results, and postpone the sketch of the derivation to Sec. 13. 31 

In the intermediate strength tunneling regime we obtained the CGF in the following form 1 46 1 



Hx) « |^m[l + T f H/ L H{l^^ (28) 

where fr(uj) = /(w— fi T ) is the Fermi distribution function for electrons in the lead r. This result looks simi- 
lar to the Levitov-Lesovik formula for noninteracting systems 1 2 1, but the effective transmission probability 
T f {uj) accounts for the strong quantum fluctuations of the charge |46 53 1, 

T f {uj) = (27T) 2 , ""f^-^^-^^ coth coth (29) 



, • w - Ur \ „ , / E c \ 2nT 
2Re ib \i — ] -2ib\ — — 



- iira r (uj - fi T ) coth - ^ , (30) 



where ib is the digamma function. For a symmetric SET at T = and \oj\ <C eV, the self-energy becomes 
Si=l r ^f^) ~ a o \n{2Ec/eV) u> — iT/2. From the real-part of this expression one reproduces the 
logarithmic renormalization J27> . and thus in this sense we accounted for leading logarithms. The imagi- 
nary part describes the effect of finite life-time of the charge state of the quantum dot in non-equilibrium 
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Fig. 6 (a) The zero-temperature current distribution at eV = 0.2Ec for various values of Ao ( i?L = Rr and Cl = 
Cr). (b) Plot of P at Ao = for various values of the conductance versus the current normalized to V/2(R\j + Rr); 
inset: the same distribution normalized to the average current (/). 



situations. We also reproduce the average current predicted by Schoeller and Schon 1 37 1 . Furthermore, the 
second order expansion in ao, 

jr(i) + jr( 2 ), with 

F {2) (X) = a Ao {ReSf (A )^ 1 )( X )}+^ cot (x), T cot ( X ) = t {j + (e^-l) +1 -(e-^-l)}, (31) 

is consistent with the result of Ref. 1 36 1. The first term in the expression for T^ 2 ' (\), describes the renor- 
malization 1271 in lowest order perturbation theory. Corrections of this type for the current were derived 
earlier in Ref. |38|. The cotunneling correction to the CGF, jF cot , describes a bidirectional Poissonian 
process governed by the cotunneling rates 

7 =2ira Q a du — — t— -r^r Re- 



(e±("-ML)/T - l)(l - 6 t("-mr)/t) (uj + M- A ) 2 ' 

This term dominates in the Coulomb blockade regime (|An| > \eV/2\ for a symmetric SET) and is con- 
sistent with the FCS theory of quasiparticle tunneling in the presence of many -body interaction |54|. 



3.2 NON-MARKOVIAN EFFECTS : RENORMALIZATION AND FINITE 
LIFETIME BROADENING OF CHARGE STATES 

We present now some results. Figure |6j a) shows the current (/ = eN/to) distribution for a symmetric 
SET and for several values of the Coulomb energy barrier Ao. The conductance is chosen to be very 
small. As we sweep Ao from the center of the conductance peak, Ao/eV = 0, to the threshold of the 
Coulomb blockade regime, Ao /eV = 0.5, the CGF gradually changes from the correlated Poissonian d25t 
to the uncorrected one (1241 . Simultaneously the current distribution widens. With further increase of 
Ao one enters into the Coulomb blockade regime, where the CGF smoothly crosses over to T cot and the 
non-Markovian co-tunneling processes becomes dominant. 

As the conductance increases, quantum fluctuations are enhanced. We find that for zqT <C A, [A = 
max(|zrjAo|, 27rT, |eV|/2)], the simple expression for CGF, T^, still holds provided the parameters are 
properly renormalized a — > z ao, A n — ► z A . The effect of this renormalization is illustrated in 
Fig-EJb), where the current distribution for A = is plotted. Since z decreases with increasing ao, the 
mean value of the current, i.e. the position of the peak in the current distribution, shifts to lower values. 
The renormalization effect can be absorbed if we re-plot the same data with the current (horizontal axis) 
normalized by the average current (/) rather than by V/2(R^ + Rr) [inset of Fig.|6jb)]. However, even 
after this procedure the three curves do not completely collapse to a single one. The remaining differences 
can be attributed to the non-Markovian effect of the broadening of the charge states due to their finite life 
time, which is described by the imaginary part of the self-energy, Eq. J30i . We observe the following trend: 
the probability for current much larger than the average value is suppressed and the current distribution 
shrinks with increasing ao. 
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Fig. 7 Panels (a-1) and (b-1) are the skewness and the kurtosis at eV = O.lEc and K for various conductance. 
Panels (a-2) and (b-2) are those normalized by the average value. 

Let us discuss the lifetime broadening effect quantitatively. At moderately large voltages, eV S> 7k = 
E c e~ 1/2a ° /2tt, and at T = the real part of the self-energy (|30|l is negligible and 2~Z r =L r ~ 
—iiraoeV. Then the CGF at Ao = reads 

Hx) « F {1) (X) ~ 4 9Q (e ix -l) + 2qn 2 a 2 (e*x/2_ e <x/2) + O (a 3 ). (32) 

It is evident from this formula that the higher order cumulants are suppressed with increasing ao due to 
lifetime broadening. 

Figures0a-1) and (b-1) show the skewness C3 and the kurtosis C4 as a function of Ao. A double-peak 
structure growing with increasing conductance is observed. In general, we found that higher cumulants 
of the current depend on the gate charge in a complicated way and depend strongly on the conductance. 
E.g., the kurtosis even changes its sign for large values of ao- For the generalized 'Fano factors' defined 
as C3/C1 and C3/C1 [Figs.|6ja-2) and (b-2)], we observe a suppression with increasing ao 

3 . 3 KELDYSH ACTION AND CGF IN MAJORANA REPRESENTATION 

In this section, we sketch our calculations 1461 1531 . Since within the Schwinger-Keldysh approach the 
calculation of CGF is formally equivalent to the calculation of the partition function l|9} on the closed time- 
path C (Fig. 13, we can apply the standard field theory methods. From the Hamiltonian \26\ we obtain the 
path-integral representation of the Keldy sh partition function l|9} following the standard procedure. Tracing 
out the electron degrees of freedom, we obtain the effective Keldysh action, which is the sum of two parts: 
the charging part S c h and the tunneling one St — S = S c h + St . They read 




Here c(i) and <f>(t) are Grassmann fields which correspond to Dirac fermionic operators c and Majorana 
fermionic operators </>, respectively. The latter operators are defined by the following relations: a + = &<f) 
and <7 Z = 2c^c — 1. The action St is equivalent to the tunneling action (I13> . defined earlier. 

A particle-hole Green's function, a r , describing the tunneling of an electron between an electrode r=L,R 
and the island, is expressed in the Keldysh space as a 2 x 2 matrix, 

- n ( a?{u) \ . r (u-Hr)E% fO -1 \ .... 
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where, — / (Air 2 R r ) is the dimensionless junction conductance. Here we have introduced a Lorentzian 
cutoff function to regularize the UV divergence. The terms higher than T r 4 describe the elastic cotunneling 
process. They can be neglected in the limit of fixed cto but large number of channels, ^> 1. In this 
case T T scales as oc l/v^ch an d thus the terms oc give only small corrections oc In order to 

derive the CGF, we introduce the counting field \ by means of the following rotation in the Keldysh space, 

aX(u>)=exp(iKrXTi/2)a T (u)) exp(-ire r xTi/2). (34) 

Here n is the Pauli matrix. The CGF takes the following form 



where S(x) is the effective action containing the rotated particle-hole Green's function J34> . Tracing out 
the c fields, we obtain a term of fourth order in <fi in the action, which means that the path integration 
cannot be performed exactly. Therefore we proceed in a perturbative expansion in «o and resum a certain 
class of diagrams. Namely we take into account the contributions from free Majorana Green's function 
and Dirac Green's function with the bubble insertions formed by particle-hole Green's J34I function and 
free Majorana Green's function. 

4 FCS and Coulomb interaction in diffusive conductors 

In this section we consider the FCS of low-dimensional diffusive conductors, such as quasi-one-dimensional 
disordered wires and two-dimensional disordered films. It has been appreciated more than two decades ago 
that the interplay of interaction and phase coherence effects in these systems may drastically affect its trans- 
port properties 1551 1561 . Initially, the conductance was a main object of study, but more recently the trend 
moved toward the study of essentially non-equilibrium phenomena, like the quantum shot noise. 

We consider short diffusive wires and films (with diffusion constant D) where the Thouless energy 
-&rh = D/L 2 is large compared to the applied voltage, i?xh ^S> eV. In these systems the Fano factor, i.e. 
the ratio between shot noise and current, S = 2\e\IF, takes the value F — 1/3 1 57 , 58 1. The condition for 
short conductors can be equivalently rewritten as td <C 1/eV, where td = L 2 /D is a typical diffusion 
time of electron through the system. Such short conductors are coherent and effectively zero-dimensional 
so that all effects of Coulomb interaction come from the external electromagnetic environment. It has 
been shown recently that the environment modifies the conductance, noise |26 29 1 and generally the 
FCS I301I3T1. 

Much less is known about the role of Coulomb interaction onto the FCS in the quasi- ID and 2D 
diffusive systems, when td ^s> 1/ eV . Under this condition the inelastic electron-electron scattering inside 
the conductor is important. This subject has recently attracted the attention in Ref. |59 60 61 1, where 
the so-called "hot electron" regime , was discussed. It is defined by td ^> te, with te being the energy 
relaxation time due to Coulomb interaction, and at the same time td <C T c - P h, with T _ p h being the 
electron-phonon relaxation time. These two conditions imply that the size of the conductor is larger than 
the energy relaxation length due to electron-electron interaction, but the energy relaxation from the electron 
subsystem to phonons is negligible. In this situation the electron distribution function relaxes to the local 
Fermi distribution with a position dependent electron temperature along the conductor. This changes the 
Fano factor F from 1/3 to \/3/4 |62|, an effect that was confirmed experimentally |58 1. 

The microscopic theory 1 63 1 of electron-electron interaction in low-dimensional disordered conductors 
predicts, however, in addition to te a further time scale, the dephasing time T0(See Table I). Both times are 
energy dependent and in the limit of good conductors g = G/Gq ^> 1, which we wish to consider, differ 
parametrically from each other (r^ <C te). It is usually believed [64 1 that classical phenomena described 
by the Boltzmann equation are governed only by the energy relaxation time te, while the decoherence 
time T0 affects essentially quantum-mechanical phenomena. Since the FCS is a classical quantity, in the 




(35) 
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Table 1 The electron scattering times for low-dimensional (ID and 2D) diffusive conductors, E — max{T, eV}. At 
T < tT 1 (V), we get t* = t*(V). 



d 


1/te 




l/r*, T>r-\V) 


1 

2 


E/g 


{E 2 /Du 2 ) 1 ^ 
{E/g)hxg 


(eV/T)W r^(V) 
In(eV/T) r E \V) 



sense that it is proportional to the number of conducting channels, one might naively expect that it crosses 
over between the coherent and the "hot electron" regime on the scale te ~ te. 

As we show below the time te is indeed responsible for a smooth crossover between the coherent and 
the "hot electron" limits if one is interested in the shot noise and the 3 d cumulant of charge. However, 
this is not the case for the higher order cumulants of charge transfer in the shot noise limit eV ^S> T. 
Moreover, in this limit the smooth crossover in the FCS does not exist. The Coulomb interaction drastically 
enhances the probability of current fluctuations for short conductors 1/eV <C td -C te. We coined for 
this range of parameters the term "incoherent cold electrons" |65 1. In what follows we will show that the 
tail of the current distribution for such electrons is exponential, P(I) ~ sxp(— j\I\to/e). The fluctuations 
are strongest for low temperatures, T <C Erh> an d me Y reach the maximum on the scale td ~ T,p(V). 
In this case 7 ~ g^ 1 ^ 2 for ID wire and 7 ~ (ln g/g) 1 / 2 for 2D film. The FCS of this type can be 
understood as the statistics of a photocurrent which is generated by electron-hole pairs excited by classical 
low-frequency fluctuations of the electromagnetic field. It is remarkable that the time scale of optimal 
current fluctuations transforms to the scale t v (T), known as a decoherence time in the theory of weak 
localization |63|, provided one identifies eV with T. Therefore, in strongly non-equilibrium situation the 
time rather than te governs the crossover in the FCS between the coherent and the "hot electron" limits. 

4. 1 MODEL AND EFFECTIVE ACTION. 

We consider a quasi-one-dimensional (ID) diffusive wire of length L and a quasi- two-dimensional (2D) 
film of size L x L, with density of states v& per spin, diffusion coefficient D and large dimensionless con- 
ductance g = AirvdDL 2 ~ d 3> 1, They are attached to two reservoirs with negligible external impedance 
which are kept at voltages ±V/2. The current flows along the z direction and we concentrate on the 
incoherent regime, maxje^, T} 3> -EVh- 

To evaluate the CGF we have used the Keldysh technique and employed the low-energy field theory of 
the diffusive transport 1 66 1 which leads to the action 

S[ X ,Q,A} = J d d rTr (v<2 - i[A, Q]^ - 2m> d (d t Q) J dt J ^(a^-A 2 ) . 

(36) 

Here A = diag(Ai(t, r), Aa(t, r)) is the 2x2 matrix in Keldysh space, where A12 stand for fluctu- 
ating vector potentials in the conductor. We assume that curl A = 0, thus neglecting relativistic effects. 
The matrix Q(r, t\, t 2 ) accounts for diffusive motion of electrons and obeys the semi-classical constraint 
Q(r) o Q(r) = S(t\ — £2). Boundary conditions are imposed on the field Q in the left (L) and right (R) 
reservoirs [67], Q\ r=R = Gr and Q\ T=L = G L (x) = e l ^/ 2 G L e- lxf3 / 2 . Here G LM are the Keldysh 
Green's functions in the leads. 

With the action d36i the CGF should be evaluated as a path integral over all possible realizations A12 
and Q. In general this is a complicated task. However, in the limit 1/g 1 the problem simplifies. 
We employ the parameterization Q = e lW Ge~ lW , WG + GW — 0. Here the field W accounts for the 
rapid fluctuations of Q with typical frequencies u> ~ eV and momenta q ~ WeV/D, while G(e, r) is the 
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slow stationary Usadel Green's function varying in space on the scale ~ L. As a first step we integrate 
out the field W in the Gaussian approximation to obtain the non-linear action S(\, G, A) of the screened 
electromagnetic fluctuations in the media. We keep only quadratic terms in S, what is equivalent to the 
random phase approximation (RPA). As second step one can integrate out the photon field A and reduce 
the problem to an effective action S c s[x, G}. Then the saddle point approximation, 8S c s[x, G]/SG = 
yields the kinetic equation for G(e, r). This program is very similar to that pursued in Ref. |66|. 

In the universal limit of a short screening radius, r _1 = {Ane 2 v^) 1 ^' 2 ^> y^eV/D, we get the following 
result 



Ses[x,G] = —gL 



2-d 



d d I 



de 
2tt 



VG e (r) 



to 
' 2 



r doj d q 



(2tt 



d+l 



hi 



-(L>q 2 ) 2 -uj 2 



(37) 



where T> u is a 2 x 2 matrix operator in Keldysh space corresponding to the non-equilibrium diffuson 
propagator, 

-l 



^(r,q) 



af3 



+ (i/4) / deTr( 7 V-7 Q G e+ W2W7 P G e -W2W 



with 7° = 1, 7 1 = f 3 . The first term in 5 e g is due to the diffusive motion of free electrons, while the 
second describes the real inelastic electron-electron collisions with energy transfer uj < max{eV, T}. 

Minimizing the action S e s under the constraint G(e, r) 2 = 1 one obtains a non-linear matrix kinetic 
equation for G(e, r). It has a structure of the stationary Usadel equation 



5v(G £ (r)VG e (r)) = [z e (r),G e (r) 



with the extra matrix collision integral X e in the r.h.s 
i s—^. f duid d 



d+l 



^2^(r,q) [ 7 «G e _ w (r) 7 ^ +^G e+ul (r)j a 



(38) 



(39) 



This kinetic equation should be supplemented by the %-dependent boundary conditions at the interfaces 
with the leads, as described after Eq. i36\ . Since we consider the FCS at low frequencies, SI <C £?Th> 
there is no time-dependent term in our kinetic equation similar to that of the usual time-dependent Usadel 
equation. In this limit the collision integral (13 9> guarantees the current conservation, divj = 0, where 
j a J de Tr ^f3G c (r)VG e (r)J . The resulting CGF, !F(x), can be found by evaluating the action d37l and 
solving the kinetic equation G e (x, r). In the absence of the field x our matrix kinetic equation reduces to 



the standard kinetic equation with a singular kernel K(u>) 



d/2-2 



in the collision integral |63 64 66 1 



To derive the action d37i we have used a local approximation, i.e. we neglected gradient corrections 
proportional to (VG ~ *C VW '. In this way we incorporate only classical effects of interaction 

into the FCS. The gradient terms would be responsible for quantum corrections to the CGF, coming from 
frequencies uj ^> max{T, eV}. They are small in the parameter 1/g and are beyond the scope of this 
article. 

So far our consideration was rather general. In the following we restrict the analysis to the most inter- 
esting shot-noise limit, eV 3> T. Then the further particular solution of kinetic equation strongly depends 
on the relative magnitude of the diffusion time td compared to the voltage dependent energy relaxation 
time te(V). In the range of sufficiently high voltages, so that 775 3> te(V), the system is driven into the 
"hot electron" regime. In this limit the electron distribution function of electrons has the from of a local 
Fermi distribution with position dependent temperature T(r) set by the applied voltage eV and differing 
from the temperature T in the leads. For smaller voltages, so that 1 /eV <C to <C te (V), the electrons are 
described by a strongly non-equilibrium two step distribution function, which results from the weighted 
average of the Fermi distribution functions in the left and right leads. We thus call this situation the regime 
of "cold electrons". The behavior of the FCS is essentially different in these two regimes and we consider 
them separately in the following subsections. 
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4.2 "COLD ELECTRON" REGIME. 

The cold electron regime is defined by relation i?Th 3> 1/te. Under this condition the collision term in 
the kinetic equation is small, and one can obtain the Green's function perturbatively around the coherent 

solution obeying the Usadel equation V 2 (G e (z)V z G £ (z)) = 0. Here < z < 1 is a dimensionless 

coordinate along the current direction. This solution was found in Ref. |68 1 and at T <C eV can be written 

as 

G e ° ( X , z) = L x (z)G L (e, x) + R x (z)G R (e) (40) 
T §vah(l — z)9 x _ . , sinh z6 x „ . , /—-> 7n 

L *M= Le x ^ ) = l^f< «* = M« + 

where u — 2e lx — 1 for energies < e < eV and u = 1 otherwise. In first order in td /tb the CGF can be 
found by substituting G° in the action d37l >. The main contribution comes from frequencies T < u> < eV. 
After some algebra we obtain 

HX) = -fag/tor) J de0 2 x (e)+T c Mx)- (41) 

Here the first term is the CGF of non-interacting electrons, and -Fcoil (x) is the correction due to electron- 
electron interaction. It reads 

^ = — J dZ J (2^ ln I 1 " (£>q 2 ) 2 + a; 2 J (42) 

o w 

n( X ,z) = -4L x (z)R x (z)e^ {l-L x {z)-R x (z) - [zL x (z) + (l - z)R x (z)} (e^-1)} 

where N u = (eV/\u\ — 1), and u* = max{i?Th, T}. The Thouless energy appears in the low frequency 
cut-off u* due to finite-size effects while T takes into account the smearing of a step in the Fermi distribu- 
tion. We also note two important properties of the function II(x, z), namely (i) H(ij, z) > for imaginary 
X = ij and (ii) n(x, z) = -P 2 {z)x 2 + 0(x 2 ) at x < 1 where P 2 (z) = §z 2 (l - z 2 )(l - z + z 2 ) > 0. 

To estimate the range of validity of the result J42l > we substitute a zero order distribution function 
/o(e) = (1 — z)fp(e — eV/2) + zfp(e + eV/2) into the collision integral, /f(e) being equilibrium Fermi 
distribution. Then one estimates the 1 st order correction to be 

, , L 2 /D [ eV dio feV\ {2 - d)/2 
t e (V) J e± uj \u J 

if e± = \e ± eV/2\ eV and e± > max{i?Th, T}. By virtue of Pauli's principle this correction may not 
exceed unity, $fm < 1, which is true only for e± > e*, where the scale e* is given by 

F W)~J (eV) 2 /g 2 E Th , forD = l 

V ; ~ \ eV e^q>{-gE Th /eV} , for D = 2. y ' 

This result shows that a simple perturbation theory is valid provided e* < max{£>rh, T}. Resolving this 
inequality we obtain the condition tb < r*, where the time scale r* is presented in Table 1. The 1 st order 
perturbation theory breaks down for higher voltages, when tjj > t* . In this situation we can still obtain 
the result up to a factor of order of unity from Eq. (I42> if we use as cut-off lo* ~ e*. 

The result (I41l - I42> with cut-off lo* — max{£ , Th, T, e*} enables us to evaluate all irreducible cumulants 
Ck — — {—i) k {d k /dx k )3~{x) °f a number of electrons transfered. There is no correction to the current on 
the classical level. The interaction correction to the noise and the 3 d cumulant is small in the parameter 
td I t~e and dominated by inelastic collisions with the energy transfer li ~ eV. On the contrary, the 
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Fig. 8 Sketch of the voltage dependence of the even cumulants of the charge transfer Ci n for n > 2 for a ID diffusive 
wire (left panel) and a 2D diffusive film (right panel). The curves 1,2,3,4 correspond to temperatures _&rh ~ Ti < 
T 2 <T 3 < Ta. 



leading contribution to the higher order cumulants is due to Coulomb interaction and it is dominated by 
quasi-elastic collisions with low energy transfers cu* < u> <C eV. Up to a numerical constant the result is 

(N) / eV \ d / 2 feV\ k-i-d/2 
C 2k , 2k+1 cx ^-(^- ) (%) , for k > 2 (45) 

where (N) 3> 1 is the average number of electrons transfered and u* — max{_Exh! T, e*}. 

The voltage dependence of the n th cumulant with n > 4 at different temperatures is sketched in Fig. [8] 
Eq. ( 145 \ shows that the n + 1st cumulant of the charge transfer is parametrically enhanced versus the kth 
one by the large factor eV/ui* ^> 1. It also follows from Eq. ( 145 ) that the higher cumulants grow with 
increasing voltage for Et\i > 1/ T * an d decay for E^h < 1/t*, where the new time scale T*(eV, T) is 
parametrically smaller than te, t* -C te- (See Table I). The current fluctuations are strongest if T < 
Et\i ~ l/ 7 >(^)- In this case their maximum occurs at eV/ETh ~ g for ID and at eV/E^ ~ g/lngfor 
2D. To clarify the physical origin of this strong amplification of the current fluctuations we present below 
a heuristic interpretation of the result ( I42> by relating it to the phenomenon of photo-assisted shot noise. 

Photo-assisted shot noise has been theoretically predicted by Lesovik and Levitov 11691 . They considered 
the mesoscopic scatterer with a single transmission channel T biased by the AC voltage V(t) — Vq sin(f2t) 
(See Fig. |3j. This voltage leads to an oscillating phase ip(t) — J_ oo eV(T)d,T — <I>QCOsfii across the 
conductor with amplitude $52 = eVn/hfl. It has been shown in Ref. |69| that such a phase modulation 
results in a zero-frequency non-transport shot noise due to the excitation of electron-hole pairs in the leads. 
At low temperatures, T <C ?i$n. the noise is 

e 2 £2? 

S 2 (o; = 0) = — T(l-T) ]T |nO|j2($ n ) (46) 



where P n = J 2 ($o) is the probability to excite an electron-hole pair with the absorption of n photons, 
and J n (x) are Bessel functions. In the limit of weak phase oscillations, $0 <C 1, the noise is quadratic in 
the amplitude, S 2 = GqUVlT{1 - T)$&. 

The physical origin of the interaction correction to the FCS in diffusive conductors has very much 
in common with the generation of photo-assisted shot noise. Exploiting the path integral formulation of 
quantum mechanics one can represent an interacting electron problem by a picture where a given electron is 
moving in a fluctuating electromagnetic field Ai^r, t) created by all other electrons (with indices j = 1, 2 
referring to a forward and a backward time evolution operator). Since the main effect of the interaction 
comes from low-frequency fields with uj <C eV, the classical part A = (Ai +A 2 )/2 is of main importance. 
The field A q w leads to the excitation of electron-hole pairs, which, similar to photo-assisted shot-noise, 



Copyright line will be provided by the publisher 



20 D. A. Bagrets, Y. Utsumi, D. S. Golubev, and G. Schbn: Full Counting Statistics of Interacting Electrons 




AC2 AC4 



Fig. 9 Diagrammatic representation of interaction corrections to the shot noise and C4. Shaded blocks correspond 
to the imaginary part of diffuson, denoted as T> q (t), and thick solid lines correspond to the screened propagator of the 
electromagnetic fluctuations B(t). Each vertex tj. brings time derivative dt k - 



produce corrections to the FCS of the form 

-2(l-cf) 



i(qi-q 2 )r 



(47) 



Up to second order in \ 2 the polarization operator Il(x, z) agrees with z) in Eq. ( I42> . It is important 
that this correction is proportional to the total conductance g of the system. Thus it may become comparable 
with the non-interacting result for the CGF when the magnitude of phase fluctuation across the sample 
becomes of the oder of unity. 

In diffusive system with large conductance g ^> 1 the fluctuation of A q w are screened and can be 
considered as Gaussian with an Ohmic spectrum B(lo) = (|A q w | 2 ) = N u / (vDw), where N u = eV/\uj\ 
is a non-equilibrium distribution function of electromagnetic modes. Thus in order to obtain the interaction 
contribution to the CGF, the correction ( I47i has to be exponentiated and averaged over these fluctuations. 
Such considerations give exactly the result 1421 with polarization operator calculated to accuracy x 2 . 

It is also instructive to consider the physical picture of photo-assisted current fluctuations in the time 
domain. To illustrate the main idea we compare the interaction corrections to the shot noise C 2 and the 4 th 
cumulant C4 for a ID wire. In the time representation they read 

AC 2 = jf P 2 (z)dz^2jv' ti (q,t 1 ~t 2 )B' t2 (t 2 -t 1 )dt 1 dt 2 (48) 
AC 4 = 6^ P|(z)dz^y"^ i ( (Z ,i 1 -t 2 )p; 4 ( g ,t4-i 3 )^ 3 (t 3 -t 1 )^ 2 (t 2 -t4)^i...^4 



where 



/3(/>=P^^^, and 5{q t t)= f -Im ( — ~)< " >49) 

l/e* 



7T vDiO Jq 7T \Dq 



The corresponding Feynman diagrams are shown in Fig. [9] The correlation time of the diffuson propagator 
T> q (t) is short, given by ry ~ 1/eV, while the photon propagator B(t) is strongly non-local in time with 
long correlation time r ~ 1/e* ^> ry . Therefore the time integral in AC 2 is dominated by the short range 
|ti — t 2 \ ~ Ty only, and the correction to shot noise is small. In contrast, when evaluating the correction 
to AC4 both short, \ti — t 2 \ ~ 1*3 — ^1 ~ Ty, and long time intervals, \ti — ta\ ~ \t 2 — t^\ ~ 1/e*, 
are essential. The same structure holds for higher cumulants C 2n> 4 which are expressed by one-loop 
diagrams with n diffusons and n propagators of the electromagnetic field. We thus see that a conversion of 
the long-time electromagnetic field correlations into the current of electron hole pairs is the reason for an 
enhancement of higher order cumulants in diffusive wires and films. 
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I/<I> 

Fig. 10 The log of probability to measure the big non-equilibrium current fluctuations (eV S> T). Curve (1), 
coherent regime; curve (2), incoherent "cold electron" regime, j = 0.2; curve (3), "hot electron" regime. 

With the results (I41II42> we can also explore the current probability distribution 

P(I)= f ^exp{-r!( x )}, n(x) = -Hx)+i(Ito/e) X . (50) 

J-7T ^ 

In the long-time limit, It /e 3> 1, this integral can be evaluated within the stationary phase approxi- 
mation. For this analysis it is important that the action T(x) has two branch points, \ — i*7> where 

7 ~ (to* / eV) 1 / 2 <C 1. The points ±17 give two threshold currents, / ± = (e/to)dS/dx\ l . which 

read (7± -</))/< J) = ±7/3. 

Provided the fluctuations are small, I~ < I < I + , the saddle point x* of the function Sl(x) lies ° n 
the imaginary axis and satisfies the condition \x*\ < 7- Thus with exponentially accuracy the probability 
distribution becomes P{I) ~ exp{— Due to the smallness of the parameter tjj /te we found that 
P(I) deviates only slightly from the probability Pq(I) of current fluctuations in the non-interacting limit. 
For larger current fluctuations, / < I~ or I > I + , the potential fl(x) does not have a saddle point any 
more, and one should use the contour Co of a zero phase, Imf2(x)| c = for the asymptotic analysis 
of the integral P(I). This contour is pinned by the branch point x = ±*7> which yields exponential tails 
in the current probability distribution 

P(7) « exp{-^(± 7 ) T -ylto/e} . (51) 

The results for the probability distribution are displayed in Fig. ^3 The Coulomb interaction does not 
change the Gaussian fluctuations, but strongly affects the tails of P(I). They describe long correlated 
"trains" in the transfered charge, in agreement with our previous discussion on the enhancement of higher 
order cumulants S n >2(D+i)- 

One can also relate the statistics ( 1421 to the photocounting statistics studied by Kindermann et al. |70|. 
In that work the FCS of incoherent radiation, which is passed through a highly transmitting barrier with 
transmission coefficient T < 1, was studied. It was shown that a highly degenerate (or classical) source 
of radiation with bosonic occupation number / U >1 produces long exponential tails in the photocounting 
distribution, P(n) oc exp(— n/ fj). The tails of the distribution 15 H are of the same bosonic nature. The 
classical electromagnetic field Aq M with lu <C eV and large occupation number N u = eV/\u>\ can excite 
electron-hole pairs with probability V = ui 2 / ((Dq 2 ) 2 + lu 2 ). This probability plays the role similar to 
the transmission coefficient T. It is enhanced due to the diffusive motion of electrons and can be of order 
of unity. The polarization operator z) describes the efficiency of the conversion of electromagnetic 
radiation into a current of electron-hole pairs. 
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4.3 "HOT ELECTRON" REGIME. 

In the "hot electron" regime, £Jj>h *C 1 /te, describing the regime of high applied voltages or long enough 
samples, the collision term in the kinetic equation (1381 dominates. Thus the saddle point solution of the 
action d37t should make the collision integral vanish. To find this solution we note the that collision term 
in the action is invariant under the gauge transformation G e (r) = e~ K '^G e (r)e K ^ r \ Here K e (r) = 
if3[7(r) + /3(r)(e — 0(r))] and 7,/3 and <fi are arbitrary functions in space. In particular, this leads to 
the conservation of the current density, divj = 0, and of the energy flow, divjg = 0, where oc 

(27r)~ 1 J edeTr ^T3G e (r)VG' c (r)^ . It is well known that the physical Green function G(e,r) with a 

local Fermi distribution / e (r) = [e( c-< ^ r )'/ T ( r ) + makes the collision term in conventional kinetic 
equations vanish. Its gauge transform, G c (r), does the same for the generalized kinetic equation fl38t . 

The four unknown functions <f>, 7, T and (3 can be found from the extremum of the simplified action 
■Shot which is obtained by substituting G e (r) into the diffusive part of the action 5 c ff. For the rest of this 
discussion we restrict ourselves to a ID wire, since for the 2D geometry shown in Fig.l, all results are 
identical to ID. We write the action in the form Shot = (27r) _1 got dzShot(z), where the spatial density 
Shot (2) reads 

2 2 

Shot (*) = -T(V 7 - W) 2 + (V 7 - W)V0 - y T 3 (V/3) 2 + — (VT 2 )V/3 (52) 

Here T(z) and cf>(z) have a meaning of a local temperature and chemical potential, while (3{z) and 7(2) 
are their quantum conjugate counterparts. The action J52I has to be minimized subject to the boundary 
conditions <K*0L o ,i= ± eV / 2 - T ( z )\ z=0 ,i= T > 1 t(0) = X and 7 (1) = /3(0) = /3(1) = 0. 

The action Shot possesses 4 integrals of motion. They are the physical current J = 5 Shot /9S7j, the 
"quantum" current M = dS hot /dV(t), the energy current J E = J<f> - 2|^T 3 V/3 + ^(VT 2 ), and the 
spatial density of the action Shot (2). Performing the Legendre transform ,we can reduce the task to the 
boundary value problem for two functions T(z) and (3(z). Since it appears not to be possible to obtain an 
analytic solution of these equations f or non-vanishing %, we solved them numerically. The results for the 
probability distribution P(I) are shown in Fig.[K)] As in the previous section it can be evaluated using the 
saddle point approximation. We can see from Fig.^|that the probability of positive current fluctuations, 
AI > 0, is enhanced in the "hot electron" limit as compared to the coherent regime, while the probability 
of negative fluctuations, AI < 0, is affected to a lesser extent. The action (1521 is equivalent to the actions 
of Refs. 1 60 6 1 1 under the appropriate change of variables. A further increase of the voltage or the sample 
size will eventually bring the system into the macroscopic regime, L ^> L c -ph- The conductor in this case 
displays only Nyquist noise S = 4ksT /R, while higher order cumulants vanish and the probability of 
current fluctuations becomes Gaussian l loTl . 

5 Summary 

To summarize, we have studied the full current statistics (FCS) of charge transfer in two important ex- 
amples of the mesoscopic conductors taking into account the effects of Coulomb interaction. First, we 
derived the FCS for a single-electron transistor with Coulomb blockade effects in the vicinity of a conduc- 
tance peak. Quantum fluctuations of the charge are taken into account by a summation of a certain subclass 
of diagrams, which corresponds to the leading logarithmic approximation. In lowest order in the tunneling 
strength our results reproduce the 'orthodox' theory, while in second order they account for renormaliza- 
tion and cotunneling effects. We have shown that in non-equilibrium situations quantum fluctuations of the 
charge induce lifetime broadening for the charge states of the central island. An important consequence is 
the reduction of the probability for currents much larger than the average value. 

We further investigated the effect of Coulomb interaction onto the FCS in one- and two-dimensional 
diffusive conductors. We have found that Coulomb interaction essentially enhances the probability of rare 
current fluctuations for short conductors, 1/eV <C to <C te, with td and te being the diffusion and 
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energy relaxation times. The fluctuations are strongest at low temperatures, T <C 1/td, and they reach 
the maximum when the sample size matches the voltage dependent dephasing length due to Coulomb 
interaction. We have shown that tails of the probability distribution of the transfered charge are exponential 
and they arise from the correlated fluctuations of the current of electron-hole pairs which are excited by the 
classical low-frequency fluctuations of the electromagnetic field in the media. 

Acknowledgements This work is part of a research network of the Landesstiftung Baden- Wiirttemberg gGmbH. 
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